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ABSTRACT 


This paper describes a development of an algorithm for assessing stand 
productivity by considering the stand variables. Forest stand productivity is 
one of the crucial information that required to establish the business plan for 
unit management at the beginning of forest planning activity. The main study 
objective is to find out the most significant and accurate variable 


combination to be used for assessing the forest stand productivity, as well as 
to develop productivity estimation model based on leaf area index. The study 
found the best stand variable combination in assessing stand productivity 
were density of poles (X2), volume of commercial tree having diameter at 
breast height (dbh) 20-40 cm (X16), basal area of commercial tree of dbh 
>40 cm (X20) with Kappa Accuracy of 90.56% for classifying into 5 stand 
productivity classes. It was recognized that the examined algorithm provides 
excellent accuracy of 100% when the stand productivity was classified into 
only 3 classes. The best model for assessing the stand productivity index with 
leaf area index is y = 0.6214x - 0.9928 with R2= 0.71, where y is 
productivity index and x is leaf area index. 
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1, INTRODUCTION 

Estimating forest stand productivity with traditional method is a challenging task, since it should be 
laborious, time consuming and costly. However, implementation of terrestrial forest inventory activities for 
determining the forest productivity in the field are often constrained by the availability of the human 
resources, cost and time. Thus, the terrestrial-based method is mostly difficult to be implemented in large 
area coverage. Many studies [1-13] proven that the variation of the standing stock is affected by its biological 
and physical factors that may be verified with leaf area index (LAI), density, basal area, volume, biomass, 
silviculture system, and site quality. 

The leaf area index is defined as one side of the leaf from the total area of leaf tissue (m7) per unit of 
land surface area without considering the shape of the leaf [14, 15]. The leaf area index (LAI) is important 
for studying the crown structure and it is widely used to describe photosynthesis and transpiration of plant 
canopy surfaces [16]. LAI also emerged as a key factor in determining the impact of global climate change of 
forest ecosystem [17] and highly correlated with growth and forest products [18]. The LAI hold a significant 
role in determining the productivity of forest stands through the role of radiation interception [19]. 
In addition, LAI is one of the important variables in the functional and physiological models of plants [20], 
large scale remote sensing models [21], as well as large models ecosystem productivity models [22]. 
Now, there are several tools have been developed for estimating LAI and assess crown structure with image 
analysis through light transmittance measurements, one of which is digital hemispherical photography 
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(DHP). The use of DHP requires complex image analysis procedures and need a specific software [23-25]. 
By knowing the LAI, we can assess crown cover conditions and productivity of a forest stand. 

Based on the above, the main objective of the study is to identify the determinant factors for 
assessing the stand productivity using the most consistent and accurate combination of biophysical variables 
such as stand density, basal area and timber standing stock. The developed discriminant functions would be 
then used to assess the productivity of the forest stands. The additional study objective is to develop a 
mathematical model for assessing the productivity class of forest stands using the LAI with DHP technique. 


2. RESEARCH METHOD 
2.1. Date and Site 

The This study was conducted within the concession area of PT Prabu Alaska, Fakfak District, West 
Papua Province from January to February 2017. Geographically, the study site is located between 132°43’ 
East longitude & 133°21’ East longitude and between 2°38’ South latitude & 3°51’ South latitude. 
Furthermore, the data processing, analysis and reporting were carried out between May 2017 and February 
2018 at the Remote Sensing Laboratory, Faculty of Forestry, IPB. 


2.2. Tools, Software, Hardware and Data 

The tools used on the ground data collection were GPS, suunto, rope, measuring tape, compass, 
range finder, flagging tape, digital camera, fisheye lens. While, the image data processing hardware used was 
a set of computer unit with the following software: Hemiview 2.1, Gap Light Analyzer v.2, ArcGIS 10.1, 
Erdas Imagine 9.1, Excel, and Minitab 17. The primary data used were identification of tree species name, 
number of trees in each observation plot, stand density, basal area, standing stock volume and hemispherical 
photos. 


2.3. Research Procedure 
2.3.1 Vegetation Analysis Data 

Field sampling technique applied in this study was a purposive sampling technique to represent the 
variation density of a crown cover and stand productivity of forest area. There are also other considerations 
such as accessibility and locations of plots were in socially allowed areas. In this study, 28 clusters with the 
size of 50 m x 50 m square-shaped were chosen purposively, in which each cluster was divided into 56 
clusters element with the grid size of 25 m x 25 m. In each of cluster element (plot), there are sub-plot having 
size of 1O m x 10 m and 5 m x 5 m. On each plot, all physical variables were measured. The measured 
variable were diameter breast height (dbh), total tree height, branch tree height, number of individuals, 
tree species names and GPS coordinate points. All tree with dbh of >20 cm was measured tn all 25 m x 25 m 
size plot, for pole size in 10 m x 10 m, and sapling measured in 5 m x 5 m. 


2.3.2 Capture and Calculation of LAI Data 

DHP point retrieval was performed at observation time in cloudy conditions, ranging from 9am to 
4pm. The camera used was Canon EOS D1200 type with 12mp resolution and hemispherical/fisheye lens 
Sigma 8mm f/3.5 EX DC. The position of the camera is facing upwards (vertically) and pointing to the north 
of the compass [16, 26] using hand levelled [26] with height of 1.8 m. Exposure and other settings are 
automatic from the camera [25, 26], after that a new shot is taken. In this research the calculation of LAI 
from DHP using algorithm software Hemiview 2.1 with threshold method [16]. 


2.4. Data Analysis 
2.4.1 Estimating Indicator of Forest Stand Productivity 

The stand indicator used as forest productivity could describe the condition of the natural forest. 
There are three indicators of a sustainable productivity examined, namely: 1) abundance of tree regeneration, 
2) residual stands and 3) standing stock potential. The indicator variables were stand density, basal area and 
stand volume. The indicators and their corresponding variables reviewed in this study are summarized in 
Table | with the following details: 
1) Individual density in the field (N) 

Individual density was obtained using the following equation: 


ae : Total of individual 
Individual density (n ha!) = ———————_ 
Plot area 
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2) Basal area (BA) 
Basal area was obtained using the following equation: 


7 (m(d)?) 


Basal area (m2ha-1) = See 


3) Volume (V) 
The volume was calculated using the following formula: 


Volume (V) (m3ha-1) = BA. Hbc. F 


where: 
V = Volume of trees (m3); BA = Basal area (m2); D = Diameter of tree at breast height(m); Hbc = Clearbole 
tree height (m); a= Phi (3.14); f= the form factor (0.7) 

The variable indicators in Table 1 are classified into 3 classes (low, moderate, high) and 5 classes 
(very low, low, moderate, high, very high). Classification is using algorithm complete linkage euclidian 
distance method. After that, a discriminant analysis function of forest stands productivity index 1s performed. 


Table 1. Forest Productivity Indicator Variables 


Indicators Variables Code 
Regeneration Sapling density (n ha‘!) X) 
Pole density of the commercial species (n ha‘!) X2 
Pole density of all species (n ha‘!) X3 
Pole volume of the commercial species (m?ha‘!) X4 
Pole volume of all species (m*ha‘!) X5 
Tree density of the commercial species with dbh >20 cm (n ha‘!) X6 
Tree density of all species diameter >20 cm (n ha’!) X7 
Basal area of the commercial species with dbh >20 cm (m~ha‘!) Xg 
Basal area of all species diameter >20 cm (m*ha‘!) Xo 
Tree volume of the commercial species with dbh >20 cm (m*ha‘!) X10 
Tree volume of all species diameter >20 cm (m*ha‘!) Xu 
Residual Tree density of the commercial species with dbh 20-40 cm (n ha’!) X12 
stand Tree density of all species diameter 20-40 cm (n ha") X13 
Basal area of the commercial species with dbh 20-40 cm (m7ha‘!) Xia 
Basal area of all species diameter 20-40 cm (m?ha‘!) X15 
Tree volume of the commercial species with dbh 20-40 cm (m?ha‘!) X16 
Tree volume of all species diameter 20-40 cm (m*ha‘!) X17 
Potential Tree density of the commercial species with dbh >40 cm (n ha‘!) Xig 
stock Tree density of all species diameter >40 cm (n ha‘!) Xi9 
Basal area of the commercial species with dbh >40 cm (m7ha-1) X29 
Basal area of all species with >40 cm (m7ha‘!) X1 
Tree volume of the commercial species dbh >40 cm (m*ha‘!) Xo 
Tree volume of all species with dbh >40 cm (m*ha‘!) X3 


2.4.2 Data Evaluation 

Prior to any further analysis, the data was evaluated descriptively by considering the data normality 
and outlier of the sample data measured in the field. Basic statistical analysis was done to obtain the mean, 
standard deviation and coefficient of variance (CV) values and then used to create indicator classes for 
estimating forest stand productivity. The multicollinearity test was also done to fulfill the statistical rules on 
discriminant analysis as well as to see the relationship between independent variables. The level of closeness 
of the relationship between the independent variables was calculated using the Pearson correlation 
coefficient. 


2.4.3 Discriminant Functions Analysis 

In this study, the discriminant function of forest stand productivity was the forward method, 
where the variables are entered firstly, secondly, thirdly etc. Then, this discriminant method simultaneously 
identifies the sequence of variables entered. 


2.4.4 Classification, Accuracy Assessment and Model Selection 

The best model selection is done with accuracy test. The best key indicators and models are 
determined based on the model with the highest accuracy of all the set variables used. The accuracy 
assessment method used a confusion matrix, particularly to calculate the overall accuracy (OA) and Kappa 
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Accuracy (KA) [28]. Model selection was done by performing significance test that expresses the significant 
difference between the Kappa accuracy of each model. These test results were then used to select the best 
model. If the difference between models has a value of z>1.96, then this means that the model differs 
significantly, otherwise 1s similar. 


2.4.5 Model Development of Productivity Index with LAI 

The forest stand productivity index was constructed using combination from the selected variable 
model (as Y variable) associated with leaf area index (as variable X) using regression analysis. 
Mathematically, the equation of puncture index function of forest stand productivity is presented as follows: 


PI = f (LAD) 


where: PI = productivity index; f (LAI) = function of leaf area index 


3. RESULTS AND ANALYSIS 
3.1. Density, Basal Area, Volume 

The density variable values were derived from the sample plots measured in the field, expressed as 
the number of trees per unit area, while the basal area and volume values were derived from the measured 
tree diameter at breast height (dbh) per unit area. The average of tree density ranging from 31.43 to 123.14 
tree ha! (Xo, X7, X12, X13, Xig, X19), while the average tree basal area is about 0.05 + 0.03 m7ha™! (X9) which 
is the basal area of all species in dbh >20 cm. The average tree volume in the dryland natural forest area of 
Papua around 0.48 + 0.40 m*’ha? (X11) which is the volume of all species in dbh >20 cm (Table 2). 

The result of coefficient of variance (CV) calculation on tree density data ranged from 44.39- 
86.81%, tree basal area ranged between 53.70-89.11%, tree volume ranged between 69.91-117.51%. 
This means that tree density, basal area and volume data has a high variation. It also matches with a 
statement that the structure of forest stands in Papua has a high diversity [7, 29]. In addition, the average 
basal area and volume showed a regular pattern based on diameter class. The value of basal area and volume 
tend to follow the diameter class. The higher the diameter class has the higher value of basal area and 
volume. The highest basal area and volume value were found in the >40 cm diameter class ranging from 0.03 
+ 0.02 m*ha™! (X2;) and 0.30 + 0.34 m*ha! (X23) for volume. The lowest basal area and volume size were 
found in the 20-40 cm diameter class ranging from 0.02 + 0.01 m’ha™ (X;s5) for basal area and 0.18 + 0.12 
m°ha‘! (Xj7) for volume. Meanwhile, the average tree density decreased by its diameter class. The highest 
density of trees in the 20-40 cm diameter class was 88.86 + 45.86 trees ha™! (X13) and the lowest in the >40 
cm diameter class by 33.71 + 27.41 trees ha"! (X19). The stand structure forms an inverted J (Figure 1) 
indicating normal growth conditions in a normal natural forest but has an inverted bell shape in a planted 
forest, that forest conditions indicated that future forest stands will be safe [1-6, 29-31]. 


Table 2. Average Mean of Forest Productivity Variables in Study Site 


Variable Mean Std CV (%) Variable Mean Std CV (%) 
X| 1992.86 631.26 31.68 X13 88.86 45.86 51.61 
Xo 391.07 190.00 48.58 X14 0.02 0.01 56.99 
X3 453.57 194.44 42.87 X15 0.02 0.01 53.70 
X4 73.83 39.37 53.32 X16 0.16 0.12 76.23 
X5 88.13 40.69 46.17 X7 0.18 0.12 69.91 
Xo 110.86 54.82 49.45 Xig 31.43 27.28 86.81 
X7 123.14 54.66 44.39 Xj9 33.71 27.41 81.29 
Xg 0.05 0.03 60.47 X20 0.03 0.02 89.11 
Xo 0.05 0.03 54.27 Xo 0.03 0.02 85.23 
X10 0.45 0.41 91.56 X02 0.29 0.34 117.51 
X11 0.48 0.40 84.41 X03 0.30 0.34 113.25 
X12 78.29 42.68 54.52 


Another factor affecting the tree density, basal area and volume is the type of tree species 
[1-6, 10-12]. The tree species found in dryland natural forests of Papua were matoa (Pometia pinnata), 
ketapang (Terminalia catappa), batu (Pterygota horsfiledii), ketaran (Koordensiodendron pinnatum), 
pala hutan (Horsfieldia irya), binuang (Octomeles sumatrana), merbau (Untsia bijuga), and other types 
(Figure 2). These tree species are the dominant species in dryland tropical rain forest in Papua [5, 7, 29]. 
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3.2. Determining Indicator Variables for Estimating Forest Stand Productivity Index 

The classification and determining of forest productivity indicator in natural forest of Papua begins 
with conducting multicollinearity test by looking at the correlation between variable X. This test 1s used in 
determining indicator of forest productivity. The result of correlation matrix obtained informs to select 11 
variables from 23 initial variables that will be used as a combination variable to estimate the indicator of 
forest productivity in Papua (Table 3). The selected variables are sapling density (X1), pole density (Xz), 
tree density >20 cm (Xe), basal area >20 cm (Xs), V >20 cm (X10), tree density 20-40 cm (X12), basal area 
20-40 cm (Xj4), V 20-40 cm (X46), tree density >40 cm (Xj), basal area >40 cm (X29), and V >40 cm (X22). 
The combination of these variables then arranged to form 2, 3, 4, and 5 variable combination with 3 and 5 
classes of forest stand productivity. This classification was aimed to obtain the most accurate model with 
simpler variables, that could be applied with cheaper cost and easily implemented in the field. 


Table 3. The best 3 Combinations using One to Five Independent Variables to Classify 3 and 5 Classes of 
Forest Productivity with Complete Linkage Dendrogram and Standardize Euclidian Distance 


Model Forest stand productivity index 3 class 5 class 
code Regeneration Residual Overall Kappa Overall Kappa 
stand Standing stock accuracy accuracy accuracy accuracy 
(%) (%) (%) (%) 
M1 X2 100 100 100 82.97 
M2 X16 98.21 96.60 100 83.83 
M3 X22 100 100 96.42 86.99 
M4 Xy X16 96.42 81.87 100 96.46 
M5 Xi4 X20 96.42 93.70 100 86.25 
M6 X16 X20 100 100 92.85 83.26 
M7 Xy Xi4 X22 100 100 98.21 85.71 
M8 X2 X16 Xig 98.21 93.36 98.21 87.94 
M9 X2 X16 X20 100 100 98.21 90.56 
M10 Xy X2 X12 X22 100 100 91.07 74.14 
M11 X) X2 Xi4 Xig 98.21 94.34 96.42 93.65 
M12 X) X2 X16 X20 100 100 98.21 90.56 
M13 Xy X2 X6 Xi4 X20 98.21 92.53 100 96.84 
M14 Xy X X6 X16 X20 96.42 90.36 98.21 96.46 
M15 Xj X2 X10 X12 Xig 100 100 100 96.17 


3.3. Discriminant Function and Productivity Index Accuracy Test 

The study found that the overall accuracy obtained from the classification of forest stand 
productivity with 3 classes and 5 classes are more than 90%. It was also shown that stand productivity 1s 
strongly influenced by stand variables such as regeneration stock (sapling and pole density), residual stand 
(density, basal area and volume dbh 20-40 cm) also standing stock (density, basal area and volume dbh 
>40cm). The most consistent and accurate indicators were used as determinants of the productivity of dryland 
natural forest stands 1.e. pole density (X2), volume dbh 20-40 cm (Xi6), and basal area dbh >40 cm (X29) 
(Table 3). 
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The accuracies provided by each five variables combination in assessing stand producitivity is 
summarized in Table 3. Even with five productivity classes, the accurate productivity assessment was 
obtained, providing accuracy of more than 85%. The single variable using individual density indicator of 
pole (X2), dbh 20-40 cm (X16) and volume dbh >40 cm (X22) provide high kappa accuracy value in 
classifying 3 classes stand productivity, 1.e., 100%, 97% and 100%, respectively. 

When classifying 5 classes of stand productivity, the study show that 3 variables are consistently 
needed, namely regeneration stock, residual stand, and stand potential. This can be seen in the 3, 4 and 5 
variables combinations, based on the test of significance to the accuracy value of the combination model, 
the combination of 3, 4 and 5 variables did not differ significantly where the z-test shows the z-value is 
<1.96. Then, it preferable to select or apply the model with high accuracy and fewer variables. The selected 
combination that could be used as a forest stand productivity index is the combination no. 9 (M9) with 
perfect accuracy of 100%. The discriminant function of the 3 classes productivity for combinations of 
individual pole density variables (X2), volume dbh 20-40 cm (Xj6) and basal area dbh >40 cm (X20) are 
presented respectively in Table 4. Productivity class I has a medium X2 value, low Xi6 and X29 value. 
Productivity class II has a high X2 value, moderate X16 and X29 value. While, productivity class III has a low 
value of X2, and high Xi6 and X20 value. 


Table 4. Discriminant Function for Classifying Forest Productivity using Variables X2, X16 and X20 


Forest Productivity Class Discriminant Function 
I Y = -5.274 + 0.022 xX + 0.002 X16 + 0.081 X20 
II Y = -33.175 + 0.0556 X2 + 0.006 Xi¢ + 0.237 X29 
Ill Y = -15.392 - 0.0058 X2 + 0.016 Xj6 + 0.479 Xo 


Remarks: I Low productivity; II Medium productivity; III High productivity 


3.4 Model Relationship Productivity Index with LAI 

As discussed above, the combination no. 9 (M9) would be the best combination of variables for 
assessing the 3 classes stand productivity having 100% overall accuracy and kappa accuracy. Based on the 
best variable regression test, the correlation between forest stand productivity index and leaf area index 
follow the linear model with the equation y = 0.6214x - 0.9928 with R* = 0.71 (Figure 3). The pattern of 
relation between 3 classes of forest stand productivity index with leaf area index shows linear model, 
it means that the increase of stand productivity value is influenced by the meeting of the vertical structure of 
the stand at the observation location. The leaf area index affects the density, basal area and stand volume per 
unit area. An increase in leaf area index will affect the standing productivity class in general. Thus, the forest 
stand productivity class could be predicted by using the leaf area index instead of using the traditional stand 
variables such as basal area, stand density and regeneration. The use of LAI variable might reduce the 
volume field measurement since the ground LAI could be measured simply by using hemispherical camera. 
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Figure 3. Regression analysis result of forest productivity index with leaf area index 
a. linear; b. exponential; c. power; d. logarithmic 
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3.5. Leaf Area Index Value 

LAI in our research was obtained from hemispherical photographs by using the threshold method in 
hemiview canopy analysis software. Each one plot represented by a hemispherical photograph which is then 
analyzed into one LAI value. Thus, in this study we have 56 LAI values which grouped based on stands 
productivity class in natural forest. Once grouped, we describe a range of LAI values on productivity class. 
Our research results indicate that there are 3 classes of productivity in study site 1.e. productivity class I, 
productivity class HI, and productivity class III. LAI interval of productivity class I 1s 2.5-3.8, productivity 
class If vary between 4.0-4.6, and productivity class III between 5.3 to 6.6. This range value showed 
minimum and maximum range of LAI values forming fairly open to very dense canopy on 3 productivity 
class in study site. These results can be used to monitor vegetation conditions based on the value of LAI in 
dryland natural forest of West Papua. Other studies indicate LAI value on broadleaf forest is 4.4-8.8 [22] and 
global averages from ground measurements of tropical evergreen broadleaf forests around 4.9 [32]. 
Furthermore, another study shows that the in situ measurements for LAI in evergreen broadleaf forest in East 
Africa vary between 3.95 to 7.4 [29] and in Amazonian evergreen broadleaf forest between 3.25 and 5.1 [34]. 
While LAI value in the agroforestry system vary between 1.0-5.3 [16]. 

Several factors are contributed to LAI value variations in productivity class. The factors closely 
related to leaf area, tree population density, and tree dimension (basal area, volume). Studies [6], [26] also 
explained that the increase in plant populations will increase the value of LAI. At the leaf level, these include 
influences associated with variation in leaf pigments, leaf internal structure, saturation, and the orientation of 
the leaf relative to solar radiation [21, 25]. Afterwards, heterogeneity in tree height, diameter and the size of 
tree gaps influences reflectance and light transmittance at the multi tree scale [6, 21, 35]. 


4. CONCLUSION 

From the foregoing analysis and discussion, the study conclude that the most consistent and accurate 
indicators used as determinants of stand productivity of dryland forest stands are pole density (X2), volume 
commercial tree dbh 20-40 cm (Xji¢6), basal area commercial tree dbh >40 cm (X20). Algorithm of complete 
linkage with standardize Euclidian distance can be used to create a classification of productivity index. It is 
also concluded that the relationship between forest stand productivity index and leaf area index follow the 
linear model with the equation y = 0.6214x - 0.9928 with R* of 0.71. The study also indicates that 
productivity variables vary in West Papua with the value of coefficient of variance (CV) calculation on tree 
density data ranged between 51.61-81.29%, tree basal area ranged between 52.82-87.88%, and tree volume 
ranged from 69.91-113.60%. Beside the above, LAI value in study site is between 2.5-6.6. Based on the 
results obtained, it 1s also concluded that the use of DHP at a height of 1.8 m could be used for estimating the 
productivity based on the LAI value, particularly in natural forest. The researchers suggest development of 
other models using DHP Techniques at some height in different ecosystem types. 
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